Observing the Onset of Effective Mass 
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The response of a particle in a periodic potential to an applied force is commonly described by an 
effective mass which accounts for the detailed interaction between the particle and the surrounding 
potential. Using a Bose-Einstein condensate of *^Rb atoms initially in the ground band of an optical 
lattice, we experimentally show that the initial response of a particle to an applied force is in fact 
characterized by the bare mass. Subsequently, the particle response undergoes rapid oscillations and 
only over timescales long compared to that of the interband dynamics is the effective mass observed 
to be an appropriate description. 
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The eoncept of the effective mass is ubiquitous in solid 
state physics, allowing for a simple semiclassical treat- 
ment of the response of a particle in a solid to an external 
force. The complex interaction between the particle and 
the surrounding potential dresses the particle with an ef- 
fective mass, distinctly different from its bare mass, and 
allows for a description of the particles dynamics based 
on Newton's second law[l]: 
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where (a) is the expectation value of the acceleration of 
the particle under an applied force F, and m'^{k) is the 
effective mass for a particle with crystal momentum k 
and band index N. The effective mass is inversely related 
to the curvature of the dispersion relation, and in ID is 
given by 
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where EN{k) is the energy of the state, and h is Planck's 
constant. The modern description of electronic conduc- 
tion in solids is intimately tied to the concept of the ef- 
fective mass. 

However, a direct application of Ehrenfest's theorem 
[2] shows that, for a particle originally in one band, the 
initial acceleration due to an applied force is F/nio, where 
mo is the bare mass, and not F/m*. This is because the 
external force unavoidably leads to interaction energies 
associated with both intraband and interband dynam- 
ics, and while the intraband portion of the interaction 
alone would lead to a response described by the effec- 
tive mass, the additional interband contribution ensures 
an initial response given by the bare mass [3, 4]. Over 
time, the interband coupling results in rapid oscillations 
in the complex amplitudes of the initial and neighbouring 
bands, and an acceleration which itself oscillates around 
F/m* (see Fig. 1). In the presence of interband dephas- 
ing these oscillations die out. The steady state however 



contains small contributions from neighbouring bands, 
as imposed by the force, such that the total acceleration 
tends to F/m* after the decay of the transients [5]. We 
use the term dynamical mass to refer to the mass asso- 
ciated with this transient response of the particle, and 
effective mass dynamics to refer to its variation in time. 
See Supplementary Information for further details on the 
theoretical description. 

In typical solid state systems, the fast timescales of the 
transient oscillations and dephasing effects have thus far 
prohibited observation of the effective mass dynamics. 
Duque-Gomez and Sipe [4] have recently revisited this 
idea specifically in the context of ultracold atoms in opti- 
cal lattices, where the narrow momentum widths and in- 
herent length and time scales involved make observation 
of long-range quantum coherent phenomena experimen- 
tally accessible. In particular, the ability to switch the 
lattice potential faster than the characteristic response 
time of an atom in that lattice allows for direct obser- 
vation of the momentum distribution, and has no analog 
in solid state experiments. Ultracold atoms have been 
used extensively to simulate condensed matter phenom- 
ena such as Bloch oscillations [6~8] and quantum trans- 
port [9-12]. 

In this letter we report on the first observation of the 
effective mass dynamics of an ultracold cloud of atoms 
in an optical lattice. By studying the particle's response 
to an abruptly applied force, we show that the initial re- 
sponse is characterized by the bare mass and is indifferent 
to the presence of the lattice. Subsequently, the particle 
response exhibits rapid oscillations at the bandgap fre- 
quency, and slow Bloch oscillations which are character- 
ized by the usual effective mass. 

We perform our experiment with a Bose-Einstein con- 
densate of ^''Rb prepared in a hybrid optical and mag- 
netic trap [13]. This trap is formed by the overlap of a 
single-beam optical dipole trap with a quadruople magni- 
etic trap. The resulting potential has cylindrical symme- 
try with radial and axial harmonic trapping frequencies 
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FIG. 1. Illustration of effective mass dynamics, (a) Lattice 
dispersion relation showing the time evolution of an initially 
single-band wavepacket under an abruptly applied force, (b) 
The corresponding acceleration (a) , normalized by the applied 
force F. The force results in an interaction energy with both 
intraband and interband contributions. The response accord- 
ing to the intraband contribution alone is characterized by 
the lattice effective mass m* and leads to Bloch oscillations 
(dashed black line). The interband contribution coherently 
couples the ground band to the excited bands. This modifies 
the dynamics such that the total acceleration (black solid line) 
is initially that of a free particle with bare mass mo (bare mass 
response shown by green dash-dotted line), and then rapidly 
oscillates around the usual effective mass behaviour. 



fr — 80 Hz and fz ~ 20 Hz, respectively. In this trap 
we produce nearly pure condensates of about 10^ atoms, 
providing us with the long-range spatial coherence nec- 
essary for observing matter-wave phenomena. We then 
ramp up a laser standing-wave pattern over 100 ms, adi- 
abatically loading our atoms into the ground state of a 
ID optical lattice with lattice constant d = A/2 = 532 
nm. The total potential has the form 



U ^ULCOs'^{krz)-F{t)z, 



(3) 



where Ul is the lattice depth, and F{t) the applied 
force (initially zero). The photon wavevector kr — 2tt/X 
sets characteristic momentum and energy scales, hkr and 
Er — h^k'^/2mo, respectively. We express the lattice 
depth in terms of the recoil energy via the dimensionless 
parameter s = Ul/Ej.. The lattice potential generated 
in this way is essentially defect-free, and thus avoids the 
complications due to scattering that arise in typical solid 
state systems. The peak atomic density in the lattice is 
less than STO^'^ atoms/cm^. In this regime inter-particle 
interactions have been shown to be negligible up to a mi- 
nor correction to the lattice depth [7]. We estimate the 
lattice depth based on the momentum sideband ampli- 
tude for the k ~ Bloch state. Comparing to a simu- 
lation of the Gross-Pitaevskii equation, we find that the 
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FIG. 2. Observation of effective mass dynamics, s — 9.4 
and F/mo — 11.7 /im/ms^. (a) Composite of raw absorp- 
tion images. Each row is a slice of an absorption image taken 
after variable evolution time in the lattice (4 ^s resolution) 
and 20 ms TOF, and is thus representative of the momentum 
distribution inside the lattice, (b) Example absorption im- 
age (upper), and ID profile and fit (lower). The velocity and 
amplitude of each peak are extracted and plotted in (c) and 
(d) respectively, (e) The average velocity is reconstructed by 
the sum of the velocity components weighted by their ampli- 
tudes. The solid black curve is the result of low-pass filtering 
this data, and serves as a guide to the eye. 



dynamics occur as if the lattice had a depth Ueff < Ul- 
Our results are thus compared to a single-particle analy- 
sis at lattice depth Ueff, which for our typical densities, 
represents a correction of less than 10%. For our typ- 
ical lattice depths and forces we can neglect interband 
Landau-Zcner tunnelling [14]. 

To initiate the effective mass dynamics, we abruptly 
shift the center of the magnetic trap by up to 1 mm. This 
shift is performed in 20 /xs, limited by the inductance of 
the coils used to generate the shift, and is much faster 
than the typical timescale corresponding to the bandgaps 
in our experiment (h/ A ^ 100 ^s). The shifted potential 
exerts a force F on our atoms along z that is essentially 
spatially uniform over the extent of our 20 ^m sample. 
After a variable evolution time in the lattice we abruptly 
switch off all trapping potentials, freezing out the lat- 
tice dynamics, and perform a 20 ms time-of-flight (TOF) 
expansion before imaging. The effective mass dynamics 
result in oscillations of the average acceleration (a) , and 
thus the average velocity (v) . The abrupt turn-off of the 
lattice preserves the momentum distribution, which we 
measure by imaging after TOF. 

Figure 2 shows the results of a typical experimental 
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run. The free-space momentum distribution (Fig. 2(a,b)) 
is a diffraction pattern consisting of components sepa- 
rated by the recoil momentum 2'hkr (the recoil velocity 
is Vr = hkr/mo = 4.3 ^m/ms). The amplitude and ve- 
locity of each peak are extracted from a fit of 4 equally- 
spaced, equal-width gaussians (Fig. 2(c,d)). The average 
velocity of the particles in the lattice as a function of 
time is then reconstructed from a weighted sum of these 
peaks (Fig. 2(e)). For each diffracted order, the momen- 
tum grows linearly in time while changing in amplitude 
as the particle travels across the Brillouin zone. As one 
component at — fc rises, a component at +k falls, result- 
ing in a periodic modulation in the average velocity. Our 
data show a clear oscillation in the average velocity on a 
millisecond timescale, consistent with Bloch oscillations 
[6, 7, 15, 16]. This phenomenon arises due to the long- 
range, inter- well coherence of the Bloch states and occurs 
at a frequency ujb = Fd/Ti [17, 18]. In addition to the 
Bloch oscillation we observe much faster dynamics on 
a 100 jjLB timescale, consistent with an oscillation at the 
bandgap frequency. This manifests itself as a modulation 
of the relative amplitudes of the diffracted momentum 
components, as expected for an interband phenomenon. 
Under the abruptly applied force, the initially single- 
band state will acquire amplitudes in adjacent bands over 
time. The coupling to these additional bands provides 
contributions to the average acceleration that will oscil- 
late at the respective band energy difference. For narrow 
momentum width wavepackets and times shortly after 
the force is applied, the acceleration is given by [3, 4] 
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where A„7v and p„N are the energy gap and momentum 
matrix element between Bloch states in bands n and 
respectively. At t = the contributions to (a) are in 
phase and the initial acceleration is that of a particle 
with bare mass, as expected from Ehrenfest's theorem, 
and can be seen in Eq. (4) by applying the effective mass 
sum rule [3]. For a particle initially in the ground band 
(iV = 1), the coupling is primarily to the first excited 
band (n = 2), and the effective mass dynamics in this 2- 
band case are governed entirely by the bandgap A21. As 
the particle traverses the Brillouin zone during a Bloch 
cycle, the bandgap changes, resulting in a variation in the 
amplitude and frequency of the effective mass oscillation. 

To quantitatively study the dynamics, we fit the the 
average velocity to the sum of two sinusoids, explicitly 
separating the Bloch oscillation from the effective mass 
dynamics: 

v{t) = Ad sin(w<jt + (t>d) + Ab sixi{u!Bt + (t>B), (5) 

where A is the amplitude, oj the frequency, and (j) the 
phase of the oscillation. The subscripts d and B indicate 
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FIG. 3. (a) Variation of dynamical mass frequency (upper) 
and Bloch frequency (lower) with applied force. The dynami- 
cal mass frequency, normalized by the bandgap A, is expected 
to be independent of the force. The Bloch frequency scales 
linearly with the applied force, (b) Variation in Bloch fre- 
quency (upper) and dynamical mass frequency (lower) with 
lattice depth. The Bloch frequency, normalized by the force 
divided by the bare mass F/mo, is expected to be a constant 
mod/h (blue solid line). The dynamical mass frequency is 
compared to the bandgap at fc = and k = kr (upper and 
lower solid blue lines, respectively). The dashed black curve 
is a comparison to simulated data fitted to Eq. 5 (see text). 



fitting parameters for the effective mass and Bloch oscil- 
lation, respectively. Due to the variation of the bandgap 
as the particle traverses the Brillouin zone, the fitting 
function used is not strictly correct at times comparable 
to the Bloch period. This fitting function is a compromise 
between capturing as many of the features of the dynam- 
ics as possible while still obtaining reliable fits. When 
extracting the effective mass oscillation we fit only the 
first 300 /iS of data (roughly 3 periods of the fast oscilla- 
tion) . The choice of 300 fj,s reduces the fitting error due 
to this variation in bandgap, while providing a sufficient 
number of cycles to obtain an accurate estimate of the 
effective mass dynamics near the center of the band at 

Figure 3 plots the dependence of these timescales on 
the applied force F and lattice depth s. The frequency 
of the slow oscillation is observed to scale linearly with 
the applied force (Fig. 3(a), lower), as expected for a 
Bloch oscillation. We also plot this frequency, scaled by 
the applied force, against lattice depth (Fig. 3(b), up- 
per) to show that the frequency is independent of lattice 
depth. The fast oscillation increases with lattice depth 
(Fig. 3(b); lower) and thus the bandgap, while being in- 
dependent of applied force (Fig. 3(a) upper). The fit- 
ted frequency is compared to the numerically calculated 
bandgap at fc = and k = kr, representing the range of 
frequencies the particle samples as it undergoes a com- 
plete Bloch oscillation. A more direct comparison to the 
data is made by fitting Eq. (5) to the first 300 /is of data 
generated from a Gross-Pitaevskii equation simulation. 
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FIG. 4. Effective mass and dynamical mass results. Tfie 
steady-state efTective mass at the center of tfie ground band 
(A*' = 1, fc = 0) associated witfi tfie Bloch oscillation is plot- 
ted as green squares. As the lattice depth increases, the band 
flattens, and the effective mass is expected to increase (dashed 
black line). The initial {t = 0) dynamical mass is extracted 
from the full dynamics (red circles) , and is expected to be con- 
sistent with the bare mass response m* /mo = 1, independent 
of lattice depth (solid blue line) . 



in the same way as we fit to the experimental data. 

From the average velocity fits we extract the initial re- 
sponse of the particle in the k = state to an applied 
force. The external force is applied within 20 /xs after 
a 20 IIS delay. This delay time is accounted for by the 
phase in the sinusoids of our fitting function. This phase 
is negligible for the Bloch oscillation but not for the effec- 
tive mass dynamics. The initial response is evaluated at 
the point in time = when the phase of the fast 

oscillation is zero. The effective mass theorem Eq. (1) de- 
scribes the response of a particle to a force over timescales 
long compared to the interband dynamics; thus we esti- 
mate this effective mass from the Bloch oscillation alone 
(green squares in Fig. 4). Since the band curvature de- 
creases with increasing lattice depth, the effective mass 
increases. However, the full response contains contribu- 
tions from both the "single-band" Bloch oscillation and 
the "multi-band" effective mass dynamics. The dynami- 
cal mass is estimated from the sum of these contributions, 
and at to is observed to be consistent with the bare mass, 
independent of lattice depth (red circles in Fig. 4). In 
deep lattices (s > 10), the deduced effective mass begins 
to deviate from the predictions. This is due to the growth 
of high-order diffraction peaks in the momentum distri- 
bution which lie beyond our imaging window. Neglecting 
these peaks causes us to overestimate the amplitude of 
the Bloch oscillation, and thus underestimate the effec- 



tive mass. Note that this has a minimal impact on the 
estimate of the dynamical mass since the dominant con- 
tribution in the deep lattice regime is from the effective 
mass dynamics. The error bars are given by the fit un- 
certainties, where the main contribution in the effective 
mass uncertainty comes from fitting the long timescale 
Bloch oscillation. In the presence of excited band de- 
cay and interband dephasing, the effective mass dynam- 
ics are expected to reduce to the behaviour described by 
the usual effective mass. In our system, this can occur 
due to inter-particle scattering [19, 20]. However, within 
the parameter range and timescales probed, we do not 
expect to see dephasing of the effective mass oscillation. 

The initial response at F/ttiq can be understood in- 
tuitively from a microscopic perspective. Consider, for 
example, a wavepacket with a narrow momentum width 
in the ground band of a deep lattice. This wavepacket 
will consist of a superposition of many small wavepackets, 
each localized at a potential minimum of the lattice. If 
an external force F is applied, the contribution from the 
lattice to the total force on each localized wavepacket will 
initially be zero. Thus the total force is simply the ap- 
plied force, resulting in an initial acceleration F/mo- As 
the localized wavepackets move away from the potential 
minima they begin to feel a force from the lattice, lead- 
ing to oscillations around the minima at the bandgap 
frequency. See Supplementary Information for further 
details. 

To our knowledge, the effective mass dynamics ex- 
plored in this work have never before been observed, de- 
spite its initial prediction nearly 60 years ago. Our work 
clarifies the role of the effective mass on short timescales, 
and is timely given the increasing interest in the effect of 
few-cycle laser pulses on solid targets. This burgeoning 
field has shown promise for high harmonic generation [21] 
and fast control of electronic devices [22, 23], where exist- 
ing theoretical treatments have yet to take into account 
the full dynamics of the effective mass [24] . 

The authors are grateful to A. Jofre, M. Sierckc and 
C. Ellenor for their work in developing the original Bose- 
Einstein condensation apparatus. We acknowledge fund- 
ing from NSERC and CIFAR. 



* rchang@physics.utoronto.ca 
[1] N. W. Ashcroft and N. D. Mermin, Solid State Physics 

(Saunders College, 1976). 
[2] D. Pfirsch and E. Spenke, Zeitschrift fiir Physik 137, 309 

(1954). 

[3] K. Hess and G. lafrate, Proceedings of the IEEE 76, 519 
(1988). 

[4] F. Duque-Gomez and J. E. Sipe, Phys. Rev. A 85, 053412 
(2012). 

[5] E. N. Adams and P. N. Argyres, Phys. Rev. 102, 605 
(1956). 



5 



[6] M. Ben Dahan, E. Peik, J. Reichel, Y. Castin, and C. Sa- 
lomon, Phys. Rev. Lett. 76, 4508 (1996). 
[7] O. Morsch, J. H. Miiller, M. Cristiani, D. Ciampini, and 

E. Arimondo, Phys. Rev. Lett. 87, 140402 (2001). 

[8] S. R. Wilkinson, C. F. Bhamcha, K. W. Madison, Q. Niu, 

and M. G. Raizen, Phys. Rev. Lett. 76, 4512 (1996). 
[9] J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, 

P. Lugan, D. Clement, L. Sanchez-Palencia, P. Bouyer, 

and A. Aspect, Nature 453, 891 (2008). 
[10] G. Roati, C. D'Errico, L. FaUani, M. Fattori, C. Fort, 

M. Zaccanti, G. Modugno, M. Modugno, and M. Ingus- 

cio. Nature 453, 895 (2008). 
[11] J. -P. Brantut, J. Meineke, D. Stadler, S. Krinner, and 

T. Esslinger, Science 337, 1069 (2012). 
[12] D. Stadler, S. Krinner, J. Meineke, J. -P. Brantut, and 

T. Esslinger, Nature 491, 736739 (2012). 
[13] Y.-J. Lin, A. R. Perry, R. L. Compton, I. B. Spielman, 

and J. V. Porto, Phys. Rev. A 79, 063631 (2009). 
[14] C. F. Bharucha, K. W. Madison, P. R. Morrow, S. R. 

Wilkinson, B. Sundaram, and M. G. Raizen, Phys. Rev. 

A 55, R857 (1997). 
[15] K. Leo, P. H. Bolivar, F. Brggemann, R. Schwedler, and 

K. Khler, Solid State Communications 84, 943 (1992). 
[16] T. Pertsch, P. Dannberg, W. Elflein, A. Brauer, and 

F. Lederer, Phys. Rev. Lett. 83, 4752 (1999). 



[17] F. Bloch, Zeitschrift fiir Physik 52, 555 (1929). 

[18] C. Zener, Proceedings of the Royal Society of London. 

Series A 145, 523 (1934). 
[19] R. Morandotti, U. Peschel, J. S. Aitchison, H. S. Eisen- 

berg, and Y. Silberberg, Phys. Rev. Lett. 83, 4756 

(1999). 

[20] M. Gustavsson, E. Haller, M. J. Mark, J. G. Danzl, 
G. Rojas-Kopeinig, and H.-C. Nagerl, Phys. Rev. Lett. 
100, 080404 (2008). 

[21] S. Ghimire, A. D. DiChiara, E. Sistrunk, P. Agostini, 
L. F. DiMauro, and D. A. Reis, Nature Physics 7, 138 

(2011) . 

[22] A. Schiffrin, T. Paasch-Colberg, N. Karpowicz, 
V. Apalkov, D. Gerster, S. Mhlbrandt, M. Korb- 
man, J. Reichert, M. Schultze, S. Holzner, J. V. Barth, 
R. Kienberger, R. Ernstorfer, V. S. Yakovlev, M. L 
Stockman, and F. Krausz, Nature 493, 70 (2013). 

[23] M. Schuhze, E. M. Bothschafter, A. Sommer, S. Holzner, 
W. Schweinberger, M. Fiess, M. Hofstetter, R. Kien- 
berger, V. Apalkov, V. S. Yakovlev, M. L Stockman, and 
F. Krausz, Nature 493, 75 (2013). 

[24] S. Ghimire, A. D. DiChiara, E. Sistrunk, G. Nd- 
abashimiye, U. B. Szafruga, A. Mohammad, P. Agostini, 
L. F. DiMauro, and D. A. Reis, Phys. Rev. A 85, 043836 

(2012) . 



Supplementary Information 



1 Expectations 

We consider a wavepacket initially confined to a single band, with nonzero Bloch state 
amplitudes in a region of the band characterized by an effective mass m*. If an external 
force F is then applied for t > 0, a naive argument about the subsequent motion would be 
based on thinking of the wavepacket as evolving according to an effective Hamiltonian 

Heff = xF. 

Were this valid the response would be simple, and would follow from a simple calculation. 
Identifying the velocity operator v, 
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from Ehrenfest's theorem we would have 
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{a{t)) = -(m v,H,ff m) = —, (1) 

where the acceleration 

would be characterized by the effective mass. Here the presence of the underlying periodic 
potential, which leads to an effective mass m* differing from the bare mass nio, obviously 
plays a crucial but implicit role. 

Yet a different simple argument conflicts with this. We imagine that before an external 
force F is applied we have a wavepacket with nonzero Bloch state amplitudes only in the 
lowest band and centered about k = 0. At least if the minima of the periodic potential 
{V{x) = V{x + d), where d is the lattice constant) are very deep, then this wavepacket will 
essentially be the sum of a number of small, localized wavepackets, each centered at a local 
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potential minimum. The expectation value of the force from the lattice on each such small, 
well-centered wavepacket vanishes. Consequently when the external force F is applied one 
could argue that any dynamical effect of the lattice could (at least initially) be neglected; the 
picture would be of each well-centered wavepacket subject to the external force as it would 
be were no lattice present, and thus one would expect the overall initial acceleration to be 
governed by the bare mass mo, 

(2) 
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where the superscript ° indicates the initial time when the force is applied. 

This argument could certainly be criticized as overly schematic, but its main fault is its 
implication that the conclusion (|2]) might be limited to deep lattice potentials and particular 
wavepackets; in fact, the result is more general than it might seem. The expectation value 
of the force due to the lattice on any wavepacket consisting of Bloch state components from 
only a single band vanishes, regardless of whether the lattice potential is deep or not. To 
see this, consider first a general state {ip) with coordinate representation iIj{x) = (xltp). The 
expectation value of the force due to the lattice is then given by 
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is the full lattice Hamiltonian. Restricting ourselves now to a wavepacket ip{x) consisting 
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solely of Bloch states within a given band n = N, 

i){x) = j dk CN{k)ipNk{,x) 



(5) 



where ipnkix) = {xlipnk) are the Bloch functions, a simple calculation gives that the expecta- 
tion value of the force due to the lattice on the wavepacket indeed vanishes, 



dV{x) 



dx 



P,Ho 



(6) 



from ([3]), since Ho lipnk) = Enik) \'ipnk)) where Enik) is the energy of the Bloch state of band 
n at fc, and 

(V'nfclplV'mfc') = Pnm{k)S{k - k'). (7) 

From the perspective of condensed matter physics this is no surprise. A standard example 
in that subject is a wavepacket evolving according to (jl]), where an initial wavepacket 
involving a single band would evolve as ip^x) — j- ip{x,t) = (a;|^/'(t)). 



ip{x,t) = j dk CN{k,t)^Nk{x), 

with simply CN{k,t) = Ciy{k) exp{—itEiy{k)/h). The velocity operator is given by 
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but because the momentum operator is diagonal in ([7]), in the expectation value |-0|V^(t)) 
the phase factors exp{—itE]\r{k)/h) are irrelevant, and we have the well-known result that 
(^ijj{t)\v\tjj{t)'^ is independent of time. From the perspective of Ehrenfest's theorem. 



;#)iip,-ffjw«)> = o 



where the commutator term vanishes because of (E]) and the fact that il){x,t) ([8]) is of the 
same form as ipi^x) ([5]), (%l){t)\v\'il){t)^ is constant precisely because the expectation value of the 
force due to the lattice on such a wavepacket vanishes, even as it propagates and disperses. 
Turning now to the imposition of an external force, where Ho — H with 



H = Ho-xF 



(9) 
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for t > 0, we still have 



v = ^ (10) 



and for an initial wavepacket consisting of Bloch states from a single band ^ the initial 
state \^') evolves to a state IV'(^)) more complicated than still, from Ehrenfest's 

theorem at t = we have 

p _ p 



ihnio m, 

in agreement with ([2]) because as the external force F is applied the expectation value of the 
force of the lattice on the wavepacket indeed vanishes. This is the argument of Pfirsch and 
Spenkd^. So the intuition leading to (E]) is correct, and that to (II]) is incorrect. 



2 Position operators 

Nonetheless, it would be harsh to completely dismiss the argument leading to ([1]), for it 
does capture some of the physics, and in a way that in fact can be made precise. To show 
this we allow for a more general wavepacket than considered in (15]), and here write 



or 



1^) = ^ / dk Cn{k) \^nk) , 



(12) 



(13) 



and we let the sum over n range over many bands. The integral over k is taken over one 
Brillouin zone; we choose the c„(fc) to be nonzero only in this zone and away from its edges. 
Then using (fT3ll the expectation value of the position operator is easily found to be 



ip* {x)xip{x)dx 



(14) 



dCn{k) dclik) 



dk 



dk 



Cn{k) 



dk + ^j dk ^nm{k)c*^{k)cm.{k). 
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The ^nm{k) are the Lax connection^ relating the periodic parts of the Bloch functions, 

.dUnk{x) 



dk 

where we write the Bloch functions as 



'^Umk{x)^mn{k), (15) 



V Ztt 



with Unk{x + d) = Unk{x)- It is natural to write (iMl) as 
where 

(^) intra ~ 2 ^ 
n 

+ / dk ^nn{k)cl{k)cn{k) 
can be identified as an "intraband" contribution, and 

^ dk ^nm{k)c*^{k)cm{k) 



dk 



I inter 



is an "interband" contribution; the prime indicates that the terms with n = m are to be 
excluded. We can formally write the two components as 

i^)intra = Yl / dkdk' cl{k) {'4)nk\Xi\ll)n,k') C„^{k') 
n,m 

(^) inter = / dkdk' cl{k) (^nfcl^el^mfc') Cm(A;')) 



and 



where 



{i)nk\Xi\ll)mk') = Sr, 



S{k-k')U{k)+z-^6{k-k') 



(16) 
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are taken as the matrix elements that define an intraband position operator Xi and an in- 
terband position operator Xg. The diagonal terms C,nn{k) must appear in (fT5|) if the Unk{x) 
are to be well-defined over the Brillouin zone and its extension; their appearance in f lT6|) will 
guarantee that terms of physical significance are independent with respect to global (i.e., 
position independent) but fc-dependent phase factors that may be chosen to multiply the 
V'nfc(a^) [2]. In higher dimensions the ^nn{k) become vector quantities and their curl is the 
Berry curvature, leading to anomalous velocity effect^, for example. In ID they are more 
benign. This all harks back to the very early work by Adams and BlountP; for a more re- 
cent presentation in this spirit and references to earlier work see Aversa and Sip^. It is 
useful to define Xnmik) = (1 — 5„m)^nm(^), so that Xnmik) = if n = m; the commutator 
[Ho,X] = hP/inio then leads to 

for m ^ where Enm{k) = En{k) — Em{k). 

Although Xg is associated with interband transitions and Xi with intraband transitions 
(see flTB]) ). and x = Xe + Xi, the operators Xg and Xj do not commut^and cannot be thought 
of as kinematically independent; this in itself points out that in general a simple separation 
of motion into "intraband" and "interband" components is not possible. 

Nonetheless, if we start with a wavepacket ([5]) initially confined to a single band and then 
subject to an external force, with the dynamics following from the Hamiltonian (Q, we can 
consider separately the initial effects of the "intraband component" of the interaction, —XiF, 
and the "interband component", —XeF, 

The initial acceleration due to the intraband component of the interaction is 

(«W)i°ntra = — ¥ / dkdk' c%{k)cN{k') {lpNk\[Xi, plU^Nk') , (19) 
where we have used ( ITOj) . Analogously, for the interband component we have 

(^(^))°ntcr = —T I dkdk' C%{k)cN{k') {^pNk\[Xe,PMNk') ■ (20) 

The matrix elements of the commutators [xj,p] and [xe,p] between two Bloch states are 
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found to 

{i^Nk\[Xr,p]\i^Nk') = , ,lJ {k - k'), (21) 

where m*^{k) is the usual effective mass for band A^, and 

{i^Nk\[Ae,pMNk') =lh(l- 5{k - k'), (22) 

which follows immediately from ( |2T1) and the fact that [x,p] = ih, with x = Xi + Xe- Using 
(ED) in (Uni), we find 

That is, if the wavepacket responded to only the intraband part of the interaction, it would 
respond with the effective mass (appropriately averaged over the wavepacket); this makes 
precise the intuition that led to ([H). 

For the full initial acceleration, however, we must also consider the interband part of the 
interaction. Combining the two contributions in (|T8|l . we have 

(a(t))° =( [ dk—\cN{k)A F = —. (24) 

We have recovered (fTTl) . and shown that while the intraband component of the interac- 
tion would alone lead to a response characterized by the effective mass, the result that the 
wavepacket response is characterized by the bare mass can be understood as arising from the 
combined effects of the intraband and interband components of the interaction. 



3 Modified Bloch states 

So the motion of a wavepacket subject to an external force involves in general a compli- 
cated interplay of interband and intraband motion. How then to push beyond the initial 
value of acceleration identified above, and identify at least approximately the full \ip{t)) equal 
to ([5]) at t = and evolving henceforth according to ([9])? 

Two of us earliei'^ addressed this problem employing the modified Bloch states introduced 
years ago by Adams and Wannier for situations where Landau-Zener tunneling can be ne- 
glectecP'^. These states, which we denote by |0nfc), are labeled by a band index n and crystal 
momentum k as are the Bloch states, but they are not eigenstates of the Ho] rather, they 
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are built in such a way that a wavepacket initially of the form 



dkhNik) \(j)Nk) 

evolves in time, according to the full Hamiltonian ([9]), as 

= / dkbj^{k,t) 



so that for t > the wavepacket is still formed by modified Bloch states with the same initial 
band index N. In other words, the wavepacket remains in the new "modified band" in the 
presence of an external force. The amplitude bi\i{k,t) is 

bN{k,t) ^ 6^(K)e-*^^("'*) 

to first order in the force, where we introduced a wavevector moving through the Brillouin 
zone, K = k — Ft/h, and a phase associated with the energy of band N renormalized by the 
diagonal part of the Lax connection, 

1 r' 

7Ar(/t,t) = - / {EN{ti + Ft'/h)-FiNN{t^ + Ft'/h))dt'. 

The modified Bloch states can be constructed from the usual Bloch states by a unitary 
transformatiorP. To first order in the force we have 

^ En'n{k) 

which clearly shows that the wavepacket is formed by a superposition of Bloch states 

where the largest contribution comes from the original band (associated with \ipNk))- The 
amplitudes in the other bands are essential to keep the wavepacket in the modified band 
(associated with |07Vfc)) for t > 0. Furthermore, a wavepacket of this form responds at all 
times with the effective mass m*j^{k), 

^Am^im = ( [ dk^-jribNikM 



dt ^ "^jv(^) 

and performs Bloch oscillations as it moves through the Brillouin zond^. In this sense the 
Bloch oscillations can be regarded as a single band effect with respect to the modified bands 
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associated with the modified Bloch states. 

Clearly, the wavepacket does not consist of states from only one of the original 

bands, and hence, were the force turned on at t = 0, its response at that time would not be 
characterized by the bare mass as shown before for a wavepacket of the form ([5]) at t = 0. How 
can we then describe the dynamics of a wavepacket initially consisting of Bloch states from 
only one of the original bands? Since the dynamics are easily described using the modified 
Bloch states, we construct our wavepacket {ipit)) in terms of them so that ip{x,t = 0) equals 
the ip{x) of This requires the inclusion of modified Bloch states with n ^ N, 

\m) = J2 I ^fc^n(fc,t)|0nfc), (25) 

with the bn{k,t = 0) chosen to satisfy the initial condition. The amplitudes in this case are 
given by 

fcjv(/«)e-*^" ifn = N 
bJk,t) ^ < / ,x , (26) 

to first order in the forc^. Since this wavepacket is initially in the original Bloch band A^, it 
will behave with the bare mass at t = 0; at later times, however, it will not simply respond 
with the usual effective mass because of the presence of additional amplitudes bn{k,t) for 
n ^ N. The expectation value of the acceleration will oscillate around the effective mass 
behaviour. Details have been presented earliei'^. 



4 Dynamics of the effective mass 

The deviations from the usual effective mass behaviour of the wavepacket (125|) for t > 
can be explained in terms of the intraband and interband contributions of the interaction 
with the force, as done in section [2] for the initial response. To illustrate this point we move 
back to the representation in terms of the usual Bloch states. From the amplitudes (126|1 we 
concludd^ that in this basis \ip{t)) can be written as ( IT2l) with Cn{k) replaced by 

cjk,t) = _ - . (27) 



For times t > we need to include the contributions of both commutators 



and 



[p, x] when applying Ehrenfest's theorem to calculate the expectation value of the acceleration 
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(see expression f lTT]) ). Accordingly, we write 
dkdk' 



I I di^df^ / / r 1 \ \ 

(a(t)) =J2 -^<('^^t)^rn{k',t) (^(^iJnk [p, Ho\ i^mk') + (^nfc | [x, p] |^^fc/ ) j . (28) 

n,m " 

As discussed in section [U the first commutator in this expression is nonzero only for 

n m, 



nk 



tprnk') = Enm{k)Pnm{k)5{k - k') 



The second commutator in fl28|) is also proportional to 6{k — k') and its diagonal terms (with 
respect to the band indices) are related to the effective mass, as shown in (12T!) and (l22l) . The 
commutation relation [x,p] = ih makes the off-diagonal components vanish, 

{ipnk\[x,p]\iJmk') = 0, for n^m. 

Using these results for the commutators in (l28l) and keeping terms to first order in the 
force, we conclude that the expectation value of the acceleration has two contributions 

m) = m)N + Ht)).o^- (29) 

The first one, given by 

A'', intra N, inter ' 

(30) 

is associated with the population of band through an intraband term, 

1 . . .,2 



and an interband term, 



{«(«)>«, ={ I dk 1 F, (31) 



(aW>.,„„=(/*(;^-;^)|c.WljF. (32) 

Recall that the population of band A^ is simply |cAr(/c,t)p = |cAr(fi;)p according to f E7|) : thus, 
the two contributions in (13T|) and ( 132|) match (a(t))°ntra ('^(^))inter at t = 0. 

Only at this initial time the acceleration is given exclusively by (a(t))^. At later times 
we need to consider the second term in fl2^ . which results from the coherence between band 
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N and the other bands. It can be written as 



(«W)coh= 5^ / -^c%{k,t)cn{k,t)EN^{k)p^^{k) + C.C., (33) 

where C.C. denotes the complex conjugate of the first term. If we ignored {a(t))^^^, assuming 
that {a{t)) ~ {a{t))jy is a reasonable approximation, we would conclude that the wavepacket 
would respond with the bare mass at all times, as occurs at t = 0: although the intraband 
term ( l3Ti) alone predicts a response with the usual effective mass, the interband term ( l32i) 
would "mask" this response and would result in a response with the bare mass at all times. 
However, once the amplitudes in bands n ^ N are included through fl33|) . the effective mass 
behaviour is "revealed" as (a(t)) oscillates around the usual effective mass according tcP 

{a{t)) ^ [ dk\c^{^)\'(^— + —,Y. f^^i%2 ^^ Kn(fcW(«:)e^^-(^'*)] Y (34) 

Here Re[-] denotes the real part of the argument and '~f]\fn{K,,t) = 7Ar(K, t) — 'yn{K,,t). 

An interpretation of this result is that the usual effective mass behaviour comes from 
the intraband part of the interaction (associated with Xj) between Bloch states belonging 
to the same band N. The initial bare mass response and the oscillatory behaviour around 
the usual effective mass come from the combined contributions of the interband part of the 
interaction (associated with Xe) between Bloch states of same band N and the effect of the 
coherence created between the amplitude in band and its neighbours. This is of course 
consistent with the picture in terms of modified Bloch states, where these states are built 
from a coherent superposition of usual Bloch states. 
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